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We analyze the loop quantization of the family of vacuum Bianchi I spacetimes, a gravitational 
system whose classical solutions describe homogeneous anisotropic cosmologies. We rigorously con- 
struct the operator that represents the Hamiltonian constraint, showing that the states of zero 
volume completely decouple from the rest of quantum states. This fact ensures that the classical 
cosmological singularity is resolved in the quantum theory. In addition, this allows us to adopt an 
equivalent quantum description in terms of a well defined densitized Hamiltonian constraint. This 
latter constraint can be regarded in a certain sense as a difference evolution equation in an internal 
time provided by one of the triad components, which is polymerically quantized. Generically, this 
evolution equation is a relation between the projection of the quantum states in three different 
sections of constant internal time. Nevertheless, around the initial singularity the equation involves 
only the two closest sections with the same orientation of the triad. This has a double effect: on the 
one hand, physical states are determined just by the data on one section, on the other hand, the 
evolution defined in this way never crosses the singularity, without the need of any special boundary 
condition. Finally, we determine the inner product and the physical Hilbert space employing group 
averaging techniques, and we specify a complete algebra of Dirac observables. This completes the 
quantization program. 

PACS numbers: 04.60.Pp,04.60.Kz,98.80.Qc 



I. INTRODUCTION 

Loop quantum cosmology (LQC) [l| is nowadays an 
active field of research, devoted to the application of the 
ideas and mathematical methods of the full theory of 
Loop Quantum Gravity (LQG) 0, i, Q to symmetry re- 
duced cosmological models. This application is useful at 
least in two respects. On the one hand, it allows us to 
learn and gain experience about issues that are still open 
in full LQG. On the other hand, in many cases the sym- 
metry reduced models already give us information about 
physical questions of interest. The first attempts to apply 
the techniques of LQG to symmetry reduced models can 
be found in Ref. Q. More recently, some homogeneous 
and isotropic models have been quantized to c omp letion 
in the LQC framework [1, 0, H, i, IliJlII, d, [ill along 
the revisited lines presented in Ref. [lj|. In particular, 
these studies provide new results about the fate of the 
classical singularities. Namely, the cosmological singu- 
larities are resolved dynamically in these models, as they 
are replaced with quantum bounces. 

In this paper we will discuss the loop quantization of a 
homogeneous but anisotropic model: the Bianchi I space- 
times in vacuo. Some preliminary analyses on the quan- 
tization of the Bianchi I model u sing Ashtekar variables 
were already developed in Refs. jlSl. [T6|. The merit for 
the first systematic attempts to construct the kinemati- 
cal Hilbert space and introduce a Hamiltonian constraint 
for the model in a loop quantization framework must be 
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granted to Bojowald pTj. However, apart from technical 
issues concerning the definition of the quantum operators 
(and the prescription adopted to incorporate the presence 
of a gap in the area spectrum of LQG), the analysis of 
Ref. [13| was not complete inasmuch as it did not pro- 
vide the physical Hilbert space, nor an algebra of Dirac 
observables. The first work that attempted to complete 
the Dirac quantization program, adapting the techniques 
presented in Ref. to quantize polymerically the grav- 
itational degrees of freedom of the Bianchi I spacetimes, 
was done by Chiou [l^. In that case, nonetheless, the 
considered homogeneous model was not in a vacuum, be- 
cause it included a massless scalar field. Here we will 
employ the same kind of techniques although we will use 
a slightly different quantization prescription [19], which 
seems more suitable to make manifest some relevant as- 
pects of the LQC approach, keep under rigorous control 
the definition of the quantum Hamiltonian constraint, 
and complete the analysis of the physical states. In par- 
ticular, an important feature of the quantization pro- 
posed here is that it immediately leads to the decoupling 
of the quantum states with zero volume [20] , so that they 
can be removed from the theory. Employing this fact, we 
will show that physical states can be described equiva- 
lently as solutions to a densitized Hamiltonian constraint. 
We will prove that the operator that represents this con- 
straint is a well defined self-adjoint operator. Further- 
more, we will explicitly construct the solutions to the 
constraint, find the physical inner product, and deter- 
mine a complete set of observables. 

Most of the homogenous models analyzed so far in 
LQC [i, 0, i, i, 0, Ilil,Jll [il, and in particular the 
mentioned work of Ref. [l8| , contain matter in the form 
of a homogeneous massless scalar field. This field, quan- 
tized in a standard (Schrodinger-like) way serves as an 
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internal time, providing a well defined notion of evolu- 
tion. Here, in the vacuum case, such a useful object is 
no longer available. Hence we will explore the problem 
of quantizing the variable that plays the role of inter- 
nal time in a polymeric way, a fact that may affect the 
conventional concept of evolution. In our case, we will 
choose as internal time one of the triad components. This 
choice will provide us with a certain notion of evolution, 
as we will discuss in Sec. VH. 

The classical vacuum Bianchi I spacetimes generically 
possess an initial cosmological singularity, to which we 
will allocate the origin of time. One of the main moti- 
vations of our analysis is to discuss what happens with 
this singularity in the loop quantum theory. We will see 
that, since the zero- volume states are totally decoupled, 
non-trivial physical states contain no contribution from 
them. In this sense the initial classical singularity is re- 
solved and disappears from the quantum theory. Fur- 
thermore, in the presented quantization there exists no 
correlation between the sectors of opposite orientations 
of the triads. Therefore, the defined evolution does not 
connect sectors corresponding to different orientations of 
the variable identified as the internal time. As a conse- 
quence, the singularity is not crossed and no additional 
branch of the universe emerges on the opposite side of it. 
The provided notion of evolution is well defined without 
the need to impose any special boundary condition to 
deal with the two orientations of the triads. 

The polymeric quantization of the system leads to a 
densitized Hamiltonian constraint that can be viewed as 
a discrete evolution equation on quantum states. This 
equation is a recurrence relation, which generically re- 
lates the projection of the state in three consecutive 
sections of constant internal time. However, when one 
reaches the origin in the quantization presented in this 
paper, it reduces to a relation on the two sections with 
the smallest possible values of the discrete time, in the 
sector of triad orientations under study. As a result, the 
physical states are in fact determined by their data on 
a single section. The vector space of these data can be 
provided with an inner product and one attains in this 
way the physical Hilbert space. 

An additional, important motivation for the analysis of 
the loop quantization of Bianchi I in vacuo comes from 
the consideration of a (much richer) family of cosmo- 
logical spacetimes that contain inhomogeneities, namely, 
the linearly polarized Gowdy model with T^-topology 
[2l| . This is an infinite dimensional model that pro- 
vides a most suitable arena for the attempt to extend 
the analysis of homogenous LQC to inhomogeneous sit- 
uations. The subfamily of homogeneous solutions within 
this model is just the classical vacuum Bianchi I space- 
times with T^-topology. From this perspective, the loop 
quantization of Bianchi I in vacuo is a preliminary step in 
order to face the quantization of the Gowdy cosmologies. 

Let us also comment that, during the writing of this 
manuscript, another work on the loop quantization of 
Bianchi I has appeared [23] which presents some similar- 



ities with our treatment. That work, carried out inde- 
pendently to ours, considers a simplified version of the 
quantization, where the corrections owing to the regular- 
ization of the inverse triad operator are not incorporated, 
and describes the time evolution using a massless scalar 
field, hke in Ref. [3. 

The main body of this paper is organized as follows. 
In Sec. II we construct the kinematical Hilbert space on 
which we define the elementary operators of the theory. 
The Hamiltonian constraint of the model is represented 
as a symmetric operator in Sec. HI, where we also show 
the decoupling of the zero- volume states. Employing this 
decoupling, we densitize the Hamiltonian constraint in 
Sec. IV, arriving to an equivalent quantum description of 
the system. The form of the densitized constraint allows 
us to decompose it in terms of one-dimensional operators, 
which are then analyzed in detail in Sec. V. The solutions 
to the constraint and the corresponding physical Hilbert 
space is determined in Sec. VI. Finally, in Sec. VII we 
discuss the results of our quantization and conclude. 



II. KINEMATICS 

As a first step towards the loop quantization of the 
vacuum Bianchi I spacetimes, we describe the model in 
terms of Ashtekar variables [ll,[2^. In principle, the def- 
inition of these variables makes use of a finite sized cell 
and a fiducial triad. Adopting a diagonal gauge, it was 
shown in Ref. [1^ that there is no physical dependence 
on the choice of fiducial triad if one defines the homoge- 
nous canonical variables for the model in a suitable way. 
For the sake of simplicity, we will then particularize the 
discussion to the choice of a diagonal Euclidean triad. 
The issue of the dependence on the coordinate cell of 
integration is more subtle [1^. Nonetheless, when the 
spatial sections of the Bianchi I cosmologies have a com- 
pact topology, there is a natural choice of coordinate cell. 
In fact, as we have mentioned in the Introduction, one 
of the motivations for our study is the potential applica- 
tion to the quantization of the homogeneous sector of the 
linearly polarized T'^-Gowdy model [2l[. Consequently, 
we will specialize our analysis to a compact three-torus 
topology, adopting the corresponding natural cell for our 
treatment, namely the T'^-cell with sides of coordinate 
length equal to 27r. In this way, one arrives to the fol- 
lowing non-trivial components of the SU{2) gravitational 
connection Af and of the densitized triad , 

A1 = |^<5f , K = j^S:, (1) 

so that {c%pj} = SnG-fdj. Here i,j = 1,2, or 3 are spa- 
tial indices, a is an internal SU (2) index, G is the Newton 
constant, and 7 is the Barbero-Immirzi parameter. 
The spacetime metric written in terms of the variables 
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Pi takes the form: 
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where {da;*} is the fiducial co-triad, with G S^, and N 
is the lapse function. 

In LQC, one adapts the techniques of LQG to symme- 
try reduced systems in order to construct the kinemat- 
ical Hilbert space [3, [13, The configuration vari- 
ables are provided by holonomies along edges oriented 
in the fiducial directions, and the momentum variables 
by triad fluxes through fiducial rectangles orthogonal to 
those directions. The holonomy along an edge of oriented 
coordinate length 27r/ii in the direction i is 



(3) 



where are the SU{2) generators proportional to the 
Pauli matrices, such that [Ti,Tj] — eij^T^' ■ The flux 
through the rectangle of coordinate area S** orthogonal 
to the direction i turns out to be 



E[S'] 



47r2 



(4) 



The configuration algebra is obtained from the sums 
of products of matrix elements of the irreducible repre- 
sentations of the holonomies, and is just the algebra of 
almost periodic functions of c* [ij, [l^] ■ This algebra is 
generated by the exponentials 
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which, using the Dirac ket notation, will be represented 
by the states The finite linear combinations of prod- 
ucts of these functions provide the analog of the space of 
cylindrical functions in LQG, and we will call it Cyl^. 
Thus, denoting |/ii,/i2,M3) = we have 



Cyls = span{|/ii,/i2,/i3)}. 



(6) 



The kinematical Hilbert space T^Kin = ®i^Kin the 
completion of the space Cyl5 with respect to the dis- 
crete inner product {pLi\^i'j) — (5^;/^; for each direction 
The states provide an orthonormal basis for 
'^kin- They are eigenstates of the operator pi associated 
to fluxes, while M^r simply shifts their label fii'. 
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Here, ^pi = V Gh is the Planck length. 

In LQG, the operator that represents the physical area 
has a discrete spectrum, with a minimum nonzero eigen- 
value equal to A = 2-\/37r7Zp[. It has been argued that, 
when one takes into account this fact, a minimum coordi- 
nate length for the edge of the holonomies is introduced 
in LQC [8|. The exact form in which such a minimum 



coordinate length must be incorporated is still under dis- 
cussion. At present, two prescriptions are considered in 
the literature [H, [11] ■ Here we will adopt the prescrip- 
tion introduced in Ref . [l^ , usually called the /I-scheme. 
One of the advantages of this prescription is that (as 
we will see) the quantum analysis of the system can be 
carried out to completion, and not just in an effective, 
non-fundamental way. Besides, this will allow us to re- 
visit some parts of the analysis presented in Ref. [l8| 
which, up to date, is the most complete discussion of the 
loop quantization of the Bianchi I model. In doing so, 
we will see that one can also learn some lessons about 
the quantization of other homogeneous systems like the 
isotropic ones. On the other hand, although the justifi- 
cation of this prescription from the viewpoint of the full 
theory of LQG is currently under investigation (26j . it 
is important to note that there are no inconsistencies or 
non-physical effects associated to it in cases with com- 
pact spatial topology like the one considered here, cases 
when a privileged coordinate cell exists [13] ■ 

According to this /I-scheme, the minimum coordinate 
length for each direction i is determined by the condition 
/i|]pi] — A, from which we arrive at the operator relation 
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Operators like vjpTf are defined in terms of pi by means 
of the associated spectral decomposition. Acting on a 
state we then get 
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Since the value Mi is state-dependent, the associated op- 
erator Np,^ generates a state-dependent minimum shift. 
To write down its action, it is convenient to relabel the 
states by reparametrizing so that the minimum shift 
becomes uniform. This is achieved by introducing, for 
each direction, a label Vi{fii) which satisfies the equation 
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whose solution is 

t'i(Mi) 



23 

357^ 



sgn(Mi)]Md^^'- 



(11) 



(12) 



With this relabeling, the basic operators have the follow- 
ing action in the domain Cyl^: 



Pi\Vt) 



3i/3Asgn(i;0]«i]'/']^^.), 
]^;. + l>. 



(13) 
(14) 



III. HAMILTONIAN CONSTRAINT 



In the considered model, only one constraint remains 
to be imposed once the diagonal gauge has been chosen: 
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the Hamiltonian constraint. In order to represent it as 
an operator, one first needs to express this constraint as 
a phase space function in terms of triads and holonomies, 
since there do not exist well defined operators corre- 
sponding to the connection components c*. This is done 
by the standard procedures of LQC, explained in detail in 
Ref. . In brief, one defines the curvature components 
employing holonomies along edges of coordinate length 
2Trp,i and regularizes the inverse of the determinant of 
the metric following Thiemann's procedure 0, i-e., ex- 
pressing it via the Poisson bracket of holonomies with the 
volume function. Applying these procedures and setting 
the lapse iV = 1, we arrive to the following form for the 
Hamiltonian constraint [l3| (integrated over the chosen 
T^-cell and still viewed as a classical function on phase 
space) 



Cbi — 



+ A2A3 




A1A3 




where 



A,: = 



Mi 



-sgn(pi)sin(^ic'), 



(15) 



(16) 



and (1/ \/\pi\) rcg is the regularized expression for 1 / vbij 
obtained via Thiemann's method. 

By the mentioned standard LQC procedures, this reg- 
ularized phase space function is represented by the oper- 
ator 



1 



1 1 



sgn(ft) 



where, for the /^-scheme that we have adopted, l//ii and 
Ap. are the operators defined in Eqs. © and p^ . re- 
spectively. Note that there is no factor ordering ambi- 
guity in the above formula, inasmuch as the operator in 
parenthesis commutes with all the others on the right- 
hand side of Eq. (fT7|) . Besides, it is easy to check that 
the states \vi) are eigenstates of the introduced operator. 
Explicitly, 
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b{vi) 



|l/3 



|l/3 



|l/3 



(18) 



In order to construct a symmetric operator Cbi repre- 
senting the Hamiltonian constraint (|15p . let us now con- 
sider the quantum counterpart of A^. From Eq. ([5]), it 



follows that all the factors in A^ depend only on pi except 
for sin(/iic'). This latter term can be represented by the 
operator 

sMjk&) = ^(A/'2f,. - AA-z^J, (19) 

which does not commute with pi . To obtain a symmetric 
operator for A^, we then proceed as follows. Since the 
operator 



(20) 



is non-negative, we can take its square root and adopt 
the symmetric factor ordering 



A,; = 



1 



2VA 



sin(^iC*)sgn(p, 
-I- sgn{pi)sm{fi,&) 



(21) 



Several comments are in order at this point. First, it is 
important to realize the presence of the factor sgn(pi) in 
the expression of A^, which was not properly taken into 
account in Ref. [ll|. As a function on phase space, it 
does not commute with sin(/2ic') under Poisson brackets, 
and hence their product as operators is not symmetric. 
Probably, its appearance had not been pointed out so far 
because, in the passage to the best studied case of (homo- 
geneous and) isotropic LQC, a simplification occurs that 
makes its role less important. Up to a constant factor, the 
purely gravitational part of the Hamiltonian constraint 
for isotropic models can be obtained from Eq. by 
identifying the three spatial directions [1, [l^ . In doing 
so, this gravitational part gets a factor of a squared sign 
that can be considered equal to the unity and ignored, 
instead of dealing with it as we have discussed for the 
anisotropic case. These two alternatives for the isotropic 
models can be understood as different choices of factor 
ordering. We have checked that a factor ordering like 
the one suggested here does not alter significantly the 
numerical results of Refs. 0, [1] (in fact, for situations 
of physical interest, the difference is below the numeri- 
cal errors). Nonetheless, even for the isotropic case our 
factor ordering may be more convenient in order to clar- 
ify certain conceptual and technical issues, like e.g. the 
decoupling of the zero-volume states or the properties 
of the solutions to the constraint near the cosmological 
singularity, as we will see later on. 

In addition, it is important to notice that the operator 
Ai annihilates the state \vi = 0), which belongs to the 
kernel of pi. Furthermore, the range of A^ does not con- 
tain the alluded state, so that its orthogonal complement 
in Hkiii is invariant under the action of A^. In particu- 
lar, this ensures that the action of the operator sgn(p,i) 
present in Aj is well defined [1^. 

The factor ordering adopted in the quantum Hamil- 
tonian constraint has some important advantages with 
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respect to that proposed in Ref. First, the con- 

straint is now a sum of products of symmetric operators, 
each defined on one of the Hilbert spaces Hl^^^ associated 
to each direction. As we will see, this facilitates the de- 
termination of observables and makes the construction of 
physical solutions straightforward, thus allowing to com- 
plete the quantization. Second, it is easy to check that 
the Hamiltonian constraint annihilates the proper sub- 
space H^jjj of states in the kernel of any of the opera- 
tors Pi. Such subspace is the completion of the subset 
of Cylg given j3y_Cyls = span{|t;i, z;2, ■ya); V1V2V3 = 0}. 

Since V = ^iy^fpTf is the volume operator, we will call 
^Kin subspace of zero-volume states. Furthermore, 
the properties of the operator commented above imply 
that the orthogonal complement of T^Kin invariant un- 
der the action of the constraint Cbi . Thus, the subspace 
of zero- volume states decouples from its complement and 
we can ignore it in the following, restricting our consid- 
erations exclusively to the subspace of nonzero-volume 
states. We will call this subspace HKm, whereas Cyl^ 
will denote the corresponding linear span of tensor prod- 
ucts of states \vi) such that none of the Wi's vanishes. 

As we will discuss in the next section, the decoupling 
of the zero- volume states allows one to describe the quan- 
tum system in a completely equivalent way in terms of a 
densitized version of the Hamiltonian constraint. More- 
over, since non-trivial physical states get no contribution 
from zero-volume states, the classical initial singularity 
disappears from the quantum theory, already at the kine- 
matical level. At least in this sense, the singularity is re- 
solved quantum mechanically, in a way similar to that 
originally suggested by Bojowald [l^ (see nonetheless 
[2^). We will consider this issue in more detail in the 
last section. 



the zero- volume states decouple is essential, because the 
kernel of the operator coincides precisely with that 

subspace. Thus, the inverse operator [l/V] is well de- 
fined once we have restricted ourselves to T^Kin- Note 
also that this densitization can be carried out exactly, 
without the need to simplify the theory by ignoring the 
quantum corrections coming from the regularized inverse 
volume operator (i.e., without replacing [l/V] by V). 

To reformulate the constraint in its densitized version, 
we introduce the following bijection in the dual Cyl^ 



(VI (VI 



(23) 



The transformed physical states are now annihilated by 
the (adjoint of the) symmetric densitized Hamiltonian 
constraint Cbi, defined as 



Its explicit form is 
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6162 + 6193 + 6263 



(24) 



(25) 



where 6^ is the symmetric operator 



6,; = 



A,: 



(26) 



This operator has the following action on the basis states 

\V^)■■ 



IV. DENSITIZED HAMILTONIAN 
CONSTRAINT 

In order to solve the quantum constraint, it proves 
convenient to recast it in a densitized form which is easier 
to analyze. One should remember that physical states 
are states annihilated by the Hamiltonian constraint Cbi 
and, in principle, they do not have to be normalizable 
in the kinematical Hilbert space TiKin- More precisely, 
we expect these states to live in a larger space, namely 

the algebraic dual Cyl^ of the dense set Cylg. We will 
denote one such state by (Vl- 

In order to densitize the quantum Hamiltonian con- 
straint in a rigorous manner, we have to invert the action 
of the operator 



(22) 



which (via Thiemann's procedure) entered in the defini- 
tion of our constraint. At this point, the observation that 



e^\v^)=-^^[f+iv^)\v^ + 2)-f-M\v^-2)], (27) 
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where 



and 



I±{vi) = g{v^±2)s±{vi)g{vi), 
s±{vi) = sgn('i;i ± 2) + sgn(ui). 



1 + i 


1 
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if «»^0, 
if V, = 0. 



(28) 
(29) 



(30) 



V. ANALYSIS OF THE CONSTRAINT 
OPERATOR 

One of the advantages of our quantization procedure 
is that, in order to study the properties of the constraint 
operator Cbi, we only need to analyze the operator 6^ on 
^kin- We carry out that analysis in this section. 
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A. Superselection 

As we see in Eq. ((27|) . 0^ is a difference operator with 
a step of two units in the label Vi. Given the definitions 
the function f+{vi) vanishes in the whole inter- 
val Vi £ [—2, 0], while f-{vi) is equal to zero for Vi € [0, 2]. 
Owing to this remarkable property, which can be traced 
back to our treatment of the factor sgn(pi) in the con- 
straint, the operator Qi does not relate states \vi) with 

> to those with Vi < 0. Therefore, Qi connects only 
states with labels Vi belonging to one of the semilattices 

/:± = {±(e,-f 2fc),/c e N}, (31) 

where 

N = N+U{0}, e^e(0,2]. (32) 

Semilattices corresponding to different values of Si or 
to different signs are not connected by the action of 9^. 
In other words, the Hilbert space Tif. , defined as the 
Cauchy completion of the set 

Cyl^, =span{|«,);i;, (33) 

with respect to the discrete inner product, is invariant 
under the action of Qi. Note that the kinematical Hilbert 
space for each direction, which is not separable, can be 
decomposed into these separable Hilbert spaces: 

nKin = ®eAnt®H-). (34) 

Because of the absence of physically relevant opera- 
tors that connect the different semilattices, the physical 
Hilbert space is then divided into superselection sectors. 
We can thus restrict our study to any specific Hilbert 
space Tit = ^iHf., with e = (81,82,63). Equivalently, 
we could construct the theory using e.g. the Hilbert space 
TiZ = (giiTi" , since the constraint is symmetric under a 
flip of sign in the label Vi, owing to the identity 

f±{-v^) = -Mv,). (35) 

Finally, we would like to emphasize that the state 
\vi = 0) is not included in any of the superselection 
sectors, since it had been removed from the kinematical 
Hilbert space. The semilattices which one might expect 
to be connected with this state under the action of the 
constraint are those corresponding to = 2, but one 
can explicitly check that Qi is indeed a completely well 
defined operator on Hf. with domain Cyl^. . Nothing 
special occurs in the case = 2 in comparison with the 
other possible superselection sectors. 

B. Spectral analysis 

In order to determine the spectral properties of the 
operator Qi it is helpful to first analyze its square. This 



squared operator is also important by itself as it rep- 
resents (up to a multiplicative constant) the gravita- 
tional part of the densitized constraint in the isotropic 
case, where the elementary variables corresponding to 
the three different spatial directions are identified. 

One can easily check that Of is a difference operator 
of constant step equal to four in the label Vi. Its ac- 
tion couples only those points Vi which lay on one of the 
semilattices 

= {±(£, + Ak), k G N}, e (0,4]. (36) 

Thus, Of leaves invariant each of the Hilbert spaces 
(4)7^^ obtained by the completion of 

(4)Cyl| = span{|t-,); e ^^^C%}. (37) 

If one now defines Of in (4)7^+ © '^^^ni_g^ (with do- 
main ('^'Cylt u('*)Cyl4„jJ, it is not difficult to check that 
its difference with respect to the operator iJ^pg A^/ (ttG) 
defined in Ref. [9] [see Eq. (37) in that reference] is just a 
symmetric, trace class operator. For the particular case 
ii — 4, we can establish the same kind of comparison with 
^APS starting with the Hilbert space (^^Tij" '^^^Ti^ 
and then including the state \vi = 0), defining e.g. a 
vanishing action of Of on it. 

Using the results obtained in Ref. [9] about the op- 
erator i?APs ^^"^ Kato's perturbation theory Q], it is 
straightforward to prove that Of is a positive, essentially 
self-adjoint operator whose essential spectrum and abso- 
lutely continuous spectrum are [0, 00) [3ll[3^. 

On the other hand, since Of leaves invariant ^^^Tif., 
its restriction to ^'^^Tif. © ^^^'H^_g. (which has been an- 
alyzed above) commutes e.g. with the projection onto 
the subspace ^'^^Tif. . As a consequence, we conclude that 

Of on the Hilbert space ^'^^Hf. (with domain ^^^Gylfj 
is essentially self-adjoint. Otherwise its deficiency index 
equation would have non-trivial solutions that would pro- 
vide also valid solutions for the case in which the operator 
is defined on the larger Hilbert space 
reaching a contradiction because we have already estab- 
lished that the operator is essentially self-adjoint in this 
latter case. Besides, for ^^"^Tif., the essential spectrum 
and the absolutely continuous spectrum must still be 
[0,00). One can show it taking into account the sym- 
metry of Of under a flip of sign in the label Vi [see Eq. 
([55)1 ] and accepting the independence of the spectrum on 
the value of e^. Moreover, numerical studies [33[ indicate 
that the whole spectrum is just absolutely continuous. 
Indeed, the spectrum of Of (with domain (^^Cyl|^) is 
non-degenerate and each of its eigenfunctions converges 
for large Vi to a non-vanishing eigenfunction of the ge- 
ometrodynamical (Wheeler-DeWitt) counterpart of the 
operator. This, together with the continuity of the spec- 
trum in geometrodynamics, suffices to conclude that the 
discrete and singular spectra are empty [3^ ]. 



7 



In a manner similar to that explained above, one can 
also relate the solutions to the deficiency index equation 
of the operators Qf and 9^ on n+ = ® ^"^H^^.^^ 

(both with domain Cylj^). In this way, one can deduce 
that the specified operator Qi is essentially self-adjoint. 
Therefore, we conclude that the constraint operator Cbi, 
given in Eq. (|25p and defined in the domain Cyli = 
^iCyl^., is in fact essentially self-adjoint. 

C. Generalized eigenstates 

Taking into account the results of the previous subsec- 
tion, we can obtain the spectral resolution of the identity 
associated to the operator Qi, e.g. on TCf., starting with 
those for the squared operator 6f on '■-^'''Hf. and . 
Remember that the spectrum of 8f on any of these two 
Hilbert spaces is absolutely continuous and equal to the 
positive real line. For all pi in the spectrum, we will call 
I ^'^^ Cp* ) the corresponding generalized eigenstates normal- 
ized to the Dirac delta t35|, where ii = Ei or 2-f e^. Thus, 
on (■*^7it we have 



(38) 



In addition, we fix the global phase of these generalized 
eigenstates by choosing ('•'^^ep'. to be a positive number 
for Vi = ii. In particular, this choice and the positivity 
of the operator 0^ ensure that i}'^^ e^^ lvi) is real for all 

Renaming pi = ojf and combining the above resolu- 
tions of the identity for the Hilbert spaces ^^"^Tif. and 
'•'^•'Ti^^g. , we obtain that, on their direct sum Hf., 



1= / dio,\etj)(e^l 



(39) 



where, for ^ 0, 



l) 



For LOi = 0, we define 




(40) 



(41) 



|e^-) = l^'^e-). 

Recalling Eq. ((27)) . one can check that the states \efj.) 

defined above are generalized eigenstates of Qi on 7i+ , 
with LUi being the corresponding eigenvalue. Whence we 
find a real, absolutely continuous spectrum. Note also 
that the projections of {efj.) on the Hilbert subspaces 
(4)-^+ ('*^7iJ"_|_^. (which are generalized eigenstates of 

the squared operator Q^ except for uji = 0, when one of 
the projection vanishes) have a relative phase of ±n/2. 
As a consequence, the phase of 



(42) 



oscillates rapidly when Vi varies in the semilattice Cf. . 

Furthermore, using Eq. (|27p . the eigenvalue equation 
associated to the operator Qi on Hf. leads to the recur- 
rence equation 



e- {2n + 2 + e.) = ^J^^-^^^^. (2n - 2 + e.) 



•UJi 



g{2n + 2 + Si) 



A g{2n + 2 + e.,)g{2n + e^) 



(43) 



Vn e , which involves three distinct values of Vi , as it 
corresponds to a second-order difference equation. How- 
ever, for n = we get a relation between the two first 
coefficients of the generalized eigenstates, 



.%/3( 



•LUi 



A g{2 + e^)g{e^) 



(44) 



Therefore the solutions to the eigenvalue problem are to- 
tally determined by a simple piece of initial data, namely, 
the projection of the generalized eigenstate for the min- 
imum allowed value of Vi, e^. (ci). Actually, a careful 
calculation shows that, Vn G N"*", 



e^_(2n + e,) = ^ 

0(n) 



n 



n 



g{2rp + Ei) 



g{2rp + 4 + E^ 



Ag{2sq + 2 + E^)g{2sq + e,) 



e!;(£.)- (45) 



Here, 0{ti) denotes the set of all possible ways to move 
from to n by jumps of one or two steps. For each 
element in 0{n), {vp} is the subset of integers followed 
by a jump of two steps, whereas {sq} is the subset of 
integers followed by a jump of only one step. 



VI. THE PHYSICAL HILBERT SPACE 

Once we have a good knowledge of the properties of 
the constraint operator, we turn to the issue of deter- 
mining the physical Hilbert space in order to complete 
the quantization. Let U be the domain of the self-adjoint 
extension of the constraint operator Cbi on the superse- 
lection sector Hi . Starting from the dense subset U, we 
can obtain the physical Hilbert space, H^^^, by applying 
the group averaging procedure [3a, [33] ■ Namely, given 
an element \(f)) in U [with corresponding wave function 
(f){v) in the ^-representation, where v — {vi,V2,V3)], one 
can "project" it onto 7i?^^ via an average over the group 
generated by the constraint Cbi: 



$(t7) = [V<j>]{v) ^ dte 



(46) 



Employing the spectral decomposition associated to 
the operators 0^ for the three spatial directions, we can 
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express the wave function 4>{v) in terms of the eigenfunc- 
tions efjXvi) introduced in Eq. 



du}cj,iu})et,\ivi)etliv2)etlM, (47) 

where uj = (wi, ^2, ij^s), (p G L^(M'^,dcD), and Vi G £f.. 
Substituting this decomposition into Eq. ((46|) . we then 
get 



— / dL05(LUlUJ2 



0^2 



(48) 



We immediately conclude from this expression that 
only products of eigenstates |e^.) with J2i(^/^i) — 
will contribute to H.^^^ . It is then useful to consider one 
of the Wi's as a function of the others. Owing to the sym- 
metry of the system under the interchange of the three 
directions, the particular choice of tOi that one makes is 
just a matter of convention. Here we select toi. From 
now on, we define it as follows: 



Wi(w2,t^3) 



LU2LO3 
UJ2 + ^3 



(49) 



With this choice, the wave function that represents the 
"projection" of the kinematical state 10) onto the physi- 
cal Hilbert space takes the form 



$(«) 



dL02duJ3 , , £ . 



2 1^2 



(50) 



where $(0^2,^3) = '^'['^1(^2, ^3), Ci;2, ^3]. We remember 
that the eigenfunctions efj.{vi) that appear in the above 
expression are given (up to a normalization factor) by 
Eq. (jiSj). 

Employing this result and the delta-normalization of 
the generalized eigenstates of Qi, one can compute the 
physical inner product between two states |$i) and |4>2) 



(*l|<&2)phy = (7'0l|<^2)Kin 
duj2dLL}3 



\UJ2 + W3I 



$t(w2,W3)«'2(w2,W3). (51) 



We have introduced the subindices "Phy" and "Kin" to 
distinguish between the inner products of the physical 
and the kinematical Hilbert spaces, and the symbol * 
denotes complex conjugation. Therefore, the physical 
Hilbert space of the considered system is the space 



Phy 



T 2 \ ^ 



1^2 +t^3| 



(52) 



of square integrable functions on with integration 
measure dLLi2duj^/\L02 + 1-^31 ■ 

An alternative way to arrive to the physical Hilbert 
space is to find the space of solutions to the constraint 



and a complete set of (real) observables, imposing then 
the condition that these observables be self-adjoint in 
order to determine a Hilbert structure on the space of 
solutions. Once we know the spectral resolution of the 
identity associated to the operators 0^, it is straightfor- 
ward to solve the constraint by adopting a formal ex- 
pansion of the states in terms of the generalized eigen- 
states |ej/. ), since Cbi has a diagonal action with this 
decomposition. Indeed, if we represent the candidate 
state by a wave function (f)'{uji,LU2,uj3), Cbi becomes 
just a polynomial constraint on LUi. Hence, the phys- 
ical solutions are described by functions of the form 
<i>'(a;2,a;3) = (j)'[uji{LU2,^3),^2,^3], where uji{lu2,^3) is 
given by Eq. (j49| . A complete set of observables is pro- 
vided e.g. by the operators O2 and O3, that multiply the 
wave function $'(u;2,'^3) by uj2 and 1^3 respectively, and 
by the derivative operators —id^^^ and —idui^. Demand- 
ing that they be self-adjoint, we arrive to the physical 
Hilbert space L'^{M.'^,du!2doJ3)- Under multiplication of 
the wave function by the factor ^y\uJ2 + i-^s], one obtains 
in fact a unitarily equivalent representation of the algebra 
of observables on the Hilbert space given in Eq. ([5^ . 



VII. CONCLUSIONS AND DISCUSSION 

In this work, we have presented a complete loop quan- 
tization of the family of homogeneous Bianchi I cosmolo- 
gies in vacuo. We have described the quantum system in 
terms of a well defined densitized Hamiltonian constraint, 
represented it by a(n essentially) self-adjoint operator, 
found the general form of its solutions, and determined 
the corresponding Hilbert space of physical states. 

A nice property of the quantization that we have put 
forward, including the chosen factor ordering, is that the 
quantum Hamiltonian constraint adopts the form of a 
sum of products of symmetric operators for each of the 
spatial directions. Under the densitization of the con- 
straint, this allows one to reduce the analysis of self- 
adjointness and the spectral analysis just to the consid- 
eration of operators for one-dimensional systems (those 
describing independently each of the spatial directions). 

Furthermore, our quantum Hamiltonian constraint an- 
nihilates the kernel of the volume operator and leaves 
invariant its orthogonal complement. As a consequence, 
the subspace of zero-volume states gets decoupled and 
one can ignore it in the study of the non-trivial physi- 
cal states. This fact is essential to attain an equivalent 
formulation of the quantum system in terms of a densi- 
tized Hamiltonian constraint. This densitized constraint 
retains all the information about the quantum behavior 
of the model; in particular, it is not necessary to intro- 
duce simplifications that disregard the regularization of 
the inverse triad operators in LQC. 

As it is common in polymeric quantizations, the phys- 
ical Hilbert space is superselected. Physical states have 
support in semilattices that have a basic step of two 
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units in the labels Vi. This differs from other previ- 
ous analysis in LQC where the constructed superselec- 
tion sectors correspond to entire lattices, with points dis- 
tributed over the whole real line instead of over a semi- 
axis 1, 0, 11 i [M [H m, El, Q [O, M Hi. As we 

have seen in Sec. V, this facilitates the spectral analysis 
of the relevant difference operators and removes remnant 
global symmetries of the system, namely those under re- 
versal of the triad orientations (which now simply relate 
different superselection sectors). In addition, and related 
with these issues, all superselection sectors have essen- 
tially the same treatment in our quantization. Nothing 
special happens for the semilattices formed by even inte- 
gers, which are those that the difference operators would 
connect with the origin had the zero-volume states not 
got decoupled. 

As we have commented, the factor ordering that we 
have adopted reduces the functional analysis of the den- 
sitized constraint operator to that of the simple opera- 
tor Qi. In addition, for the three possible directions i, 
these operators are immediately identified as Dirac ob- 
servables. It is worth commenting some of their prop- 
erties in more detail. We have seen that Qi behaves 
as a second-order difference operator in many aspects. 
Nonetheless, we know that it represents the classical 
variable c'pi, which would become a first-order differen- 
tial operator in a standard canonical quantization of the 
Wheeler-DeWitt type. Therefore, one would expect that, 
if the results of the Wheeler-DeWitt approach are to be 
recovered in a certain regime from the loop quantization, 
the generalized eigenstates have to be determined from 
data on a single section of constant Vi . Remarkably, this 
is indeed the case. 

This issue had not been pinpointed until now because 
in the homogeneous and isotropic models a subtle coin- 
cidence occurs. In that case, the gravitational part of 
the constraint is given by an operator like 8f . Such an 
operator, which represents the classical variable (c*p,i)^, 
behaves in fact as a second-order one in the Wheeler- 
DeWitt quantization. 

Another difference with respect to previous analyses 
in LOG like those presented in Refs. [1, 0, i, i, IH, 
[Til . [13 . [Tsl . 0, [1^, ll^l is that we have not introduced 
any matter content (and specifically no scalar field) as 
internal time. The variable identified with the internal 
time in our model is one of the triad components, which 
are quantized polymerically. Actually, the three diag- 
onal triad components have an equivalent role because 
the model is symmetric under the interchange of spa- 
tial directions. We have selected the direction i — 1 
to make the discussion explicit, but the choice is just a 
question of convention. Since the eigenstates of 0i are 
determined by their initial data at vi = £1, the physical 
states get also fixed by data on this initial slice. Further- 
more, the non-degeneracy of the spectrum of 0i implies 
that the wave function at any i;i-constant slice deter- 
mines the entire solution up to corrections of vanishing 
physical norm. Also the physical inner product can be 



rewritten in terms of data on that single section . One 
can then admit the viewpoint that physical states evolve 
from this initial slice to any other slice, in such a way 
that they solve the dcnsitized Hamiltonian constraint. 
As a consequence of the loop quantization adopted for 
the internal time, this concept of evolution differs from 
the usual one and, in particular, does not allow to reach 
unitarity in a straightforward, standard way ^33j] . One 
should be aware that the notion of evolution used here 
refers only to the fact that there is a deterministic re- 
lation between data on two slices of constant internal 
time, and thus its meaning is certainly limited. Also, a 
preliminary analysis of the eigenfuctions shows that the 
evolution parameter (internal time) is not monotonic in 
the cosmic time, which makes the extraction of relevant 
physical data non-trivial. 

On the other hand, the decoupling of the zero-volume 
states ensures that non-trivial physical states have no 
contribution from the slice where the curvature singular- 
ity is located in the classical theory. Physical states are 
well defined everywhere and, in this respect, the singu- 
larity is resolved quantum mechanically. This conclusion 
reinforces previous results about singularity resolution in 
LQC. It is worth emphasizing that this conclusion is in- 
dependent of our choice of internal time. Furthermore, 
since the wave functions of the physical states (in the vi 
representation) have support just in semilattices, the evo- 
lution docs not connect them with other branches of the 
universe, corresponding to a different orientation of the 
triads. The singularity is never crossed in the evolution. 
Let us also notice that this result is achieved without ap- 
pealing to any kind of boundary conditions that might 
restrict the initial data for the physical states. As far as 
one understands the statement in this sense, one can say 
that physical states "arise from nothing" in the initial 
slice vi = Si, attaining a no-boundary description. 

Apart from the points addressed in this work about the 
loop quantization of the vacuum Bianchi I model, there 
are other interesting issues that we plan to explore in a 
future research [33|. This includes a detailed compari- 
son with the Wheeler-DeWitt theory, the discussion of 
the concept of evolution, and the fate of unitarity in the 
model. We also want to carry out a numerical study of 
the physical states, analyzing in particular the occurrence 
of quantum bounces. 
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